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^ : Abstract 



The effect of collective potentials on pion spectra in ultrarelativistic heavy 



ion collisions is investigated. We find the effect of these potential to be very 
small, too small to explain the observed enhancement at low transverse mo- 



menta. 
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1. Introduction 



The enhancement at low transverse momentum found in negative particle and neutral 
pion spectra from ultra relativistic heavy-ion collisions 0, H has recently received 
considerable interest in the literature. One explanation for instance 0, [5] has interpreted 
the enhancement as a collective flow effect, but more detailed considerations concerning 
the freeze out surface result in too small an effect 0. The decay of excited baryons also 
gives rise to additional soft pions [0, || ||. However, at CERN-energies there are simply 
not enough baryons present in the central region in order to account for the observed 
enhancement ||. Furthermore the combined spectrum of thermal plus decay pions still 
flattens out at very low p t contrary to the data. 

Kataja and Ruuskanen [|1(J have shown that the measured spectrum can be fitted by 
assuming that the pions are strongly out of equilibrium. Although one may give some 
qualitative arguments [[□]] why an excess of pions should build up during the expansion, 
currently there is no quantitative understanding on how a chemical potential of the size 
needed to fit the data (/i ~ 130 MeV) would arise. 

In refs. [|12] and |13| a kinetic model with Bose-statistics in the collision integral has 
been used to study the build up of the low pt enhancement assuming different initial 
conditions for the expanding source. The authors find that the Bose-statistics in the 
collision integral indeed leads to an enhancement provided the initial pion density is 
sufficiently high. They also point out that corrections of the cross section due to Bose 
statistics are large and that they limit the effect of the Bose-phase space factors in the 
collision integral substantially. Using the corrected cross section and allowing for higher 
resonances in the initial state the authors have to assume a hadronization time as short as 
ro ~ lfm in order to account for the measured enhancement by the hadronic scattering 
processes only. 

Finally the pions may interact collectively with the surrounding hadronic medium. As 
proposed by Shuryak ]14|, H| these interactions may give rise to a strongly momentum 
dependent optical potential which is attractive for low momenta and, therefore, could lead 
to an enhanced soft component in the pion spectrum. 

It is the purpose of this article to study the effect of this latter mechanism in detail. 
In the first part we will develop the mean field potential and demonstrate how it can be 
used in a transport theoretical framework. Then we will study the effect of the mean 
fields for a static potential. Finally the expansion of the fireball will be taken in account. 
This will be done using a transport model which combines the propagation of particles in 
the mean field as well as the collisions among the particles. 
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2. Mean Field Potential 



Following Shuryak ||14|| , the most important contribution to the collective potential 



felt by a pion is the coherent scattering of two pions via the P-wave p-resonance; in the 
S-wave the contributions from different isospin channels cancel each other. We first show 
that these collective potentials can be derived from an energy functional by differentiating 
with respect to the distribution function. As a consequence energy conservation will be 
guaranteed and the collective potentials can be incorporated similar to the nuclear mean 
field in transport models of heavy ion collisions |16|, [I7| . 

For an arbitrary resonance the energy functional has the following form: 

r d 3 p d 3 p' „ , Vf , 47^ (y/i-M)r(g) 

^' J i (27r)32a;(p)(27r)32a;(p0r - Q (v^ - M? + T\q) /4 [ ) 

with q and yfs being the cm. momentum and energy and M and V are the mass and 
the width of the resonance under consideration e.g. the p. The phase-space distribution 
fi{x,p) contain an index which refers to the internal quantum numbers such as spin and 
isospin. These quantum numbers are summed over. The degeneracy factor gj j is given 
by 



(2J pl + l)(2J pl + 1) (2J p2 + l)(2/ p2 + 1) 



where the indices pi and p2 refer to the particles forming the resonance. For example in 
the case n + n — > p we would have gj t j = 1/2 

The resulting mean field potential for a particle of given internal quantum number I 
is then given by 

UAv — -a, + W f—W—Tf (x p ^J^Ml!t_ f3 ) 
Ul[P ' X) - Sfr 91 '' MP) J WMp')^ tm[ ' P) q (^-M)HP( 9 )/4 (3) 

which in the special case of 7r + 7r — > p leads to the result already obtained by Shuryak. 
Assuming isospin symmetry we obtain: 

U(p) = -*-f d3p ' fix p') AnV ~ S WL^IM (4 ) 

where f(x,p) = l/3J2i=+,-,ofi(. x >P) stands for the isospin averaged pion phase-space 
distribution. 

In the following we will truncate the explicit energy dependence of the potential U by 
evaluating (R) on the mass shell only, i.e. uj 2 = p 2 + m 2 , which is equivalent to taking 
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only the first contribution in the Dyson series for the self energy. As a consequence the 
potential depends on the 3-momentum of the particle only and, therefore, the pion wave 
function will not be modified. Hence we can treat the mesons as quasiparticles and a 
transport theoretical approach is possible^. 

For local equilibrium the distribution function can be written as 

f(x,p) = p exp(-E/T(x)) (5) 

where E = \fp 2 + m 2 denotes the energy and T(x) the local temperature. Using the 
above phase-space distribution (§) we can calculate the mean-field potential U(p) from 
eqns. (|) or (||). 

In figure [l] the resulting pion collective potential due to tt + tt — > p are shown for three 
different temperatures (T = 150, 200, 250 MeV) (the full lines). We find the potential 
becomes deeper with increasing temperature, because at the same time the density of 
pions in the heat-bath increases. The point where the potential changes sign, on the 
other hand, is essentially unaffected by the temperature. 

Since we are interested in a more general discussion of mean field effects in this article, 
in fig we show the dependence of the collective pion potential on the resonance mass. 
These potentials have been obtained for the same temperature T = 200 MeV but with 
different resonance mass M res = 500, 770, 1100 MeV. For the width we have taken a 
p-wave parameterization with a value of Tq = 150 MeV on resonance and the degeneracy 
factor was chosen to be the same as in tt + tt — > p. As we would have expected from eqn. 
(0) the potentials changes sign at a momentum close to the resonance mass. 

Since in the transport theoretical calculation described below the mean field is evalu- 
ated assuming a thermal momentum distribution, it is useful to parameterize the above 
potential in a simple form 

U(p) = V (1 - (^) 2 ) exp(-( P /a 1 ) 2 ) (^) 3 (6) 
ai T 

where Vq is the value of the potential at zero momentum while a\ corresponds to the 
momentum where the potential changes sign. The temperature dependence in this pa- 
rameterization reflects the fact that the optical potential is essentially proportional to the 
pion density. In figure |l| we compare this parameterization for three different tempera- 
tures (T = 150, 200, 250 MeV) with the potential based on tt + tt —>■ p. The parameters for 
the fit are Vq = —.2m ni a\ = Am n and To = 200 MeV. The agreement is reasonable over 
the whole range of temperatures displayed. For the sake of simplicity in the following we 

3 Note that the formalism described here may used as well for pions interacting with nucleons to study 
in-medium effects on pions in lower energy collisions. 
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Parameterset 


V [MeV] 


a x [MeV] 


T [MeV] 


Set 1 


40 


650 


200 


Set 2 


100 


650 


200 


Set 3 


100 


1300 


200 



Table 1: Different sets of parameters of effective potential (|6]). 

will only use the parameterization (g) for the mean field potential. Also the meaning of 
the parameters becomes more transparent with this choice. 

In following we will use several parametersets for the potential @ which we display 
in table [L| Set 1 corresponds to the potential obtained by Shuryak for the pure pion gas. 
Since our purpose is a general understanding of mean field effects, we have also calculated 
with other parameter sets shown in the table. Both these are much stronger than the 
realistic potential, in order to see the effects of a strong potential. 



4 



3. Static Potential 



Before we turn to the expansion of the fireball it is instructive to first study the simpler 
case of particles leaving a static potential well. Because the energy E of the particle is 
conserved while traversing the potential the energy spectra dN/ dE inside and outside are 
identical 



dN _ dN 

dE inside dE outside 

Let us assume that the particles inside the potential are distributed according to a 
boost-invariant fire tube 



dN 



dyd 2 p± 



/+oo 
drj m± cosh(^) exp(— (3mx_ coshr/) = 2m_ L K 1 (/3m_|_) 
-co 



where for simplicity we have assumed that the particles are distributed according to 
Boltzmann-statistics (f3 = Thus, in a given rapidity bin, the problem reduces to a 



two dimensional one, with E = m±_ = ym 2 + p\. 

Because of relation (|7]), the potential only affects the momentum spectra; for cylindri- 
cal symmetry we have: 



dN ldN 1 dE dN 



d 2 p p dp p dp dE 

Outside the potential we have 

dE _ p 
"dp ~ ~E~ 

thus we find 

dN 1 dN 1 dN 1 / p dN 



(9) 



(10) 



111 



inside 



d 2 p± ou t S ide E dE outside E dE inside E \dE /dp d?p 
finally using eqn. (|8|) for the momentum distribution inside the potential we find 

dN_ = Pin dN ^ 

dE outside (dE I dp)i n d 2 p± insic i e 
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where pi n is determined from 

yJpl + m* + U(j> in ) = E (13) 

and 

dE _ pi n + dU_ 
dpin \/pf„ + m 2 dpi n 



In fig |3| we have plotted the resulting spectra for the different parameterizations (table 
[l]) of the potential together with the free one. For the spectra inside the potential we have 
used a temperature of T = 1/(3 = 135 MeV, which fits the transverse momentum spectra 
from proton-proton collisions (see below). All spectra are normalized such that they have 
the same number of particles inside the potential. Therefore, in fig. |3| the integral over 
the free spectrum is larger than the integral over the ones with a potential, because in 
the latter case all particles with E < m are bound inside the well. 

We find that, as a result of the potential, the spectra become somewhat steeper at 
small momenta. The effect, however, is much too small in order to explain the data. Even 
if we increase the depth of the potential well from Vq = —40MeV (short-dashed line) to 
Vo = — lOOMeV (long-dashed line), the slope of the spectrum is not changed much (aside 
from an overall shift downwards due to the fact that the deeper potential keeps more 
particles bound inside). The difference between parameterset 2 (short-dashed line) and 
3 (long-dashed line) is even smaller. 

This behaviour actually can be understood by looking at figure |], where we have 
plotted the resulting dispersion relations E(p) for the three parameterizations. Only 
particles with energy above the rest mass E(p) > m contribute to the outside spectrum. 
There the slopes of the dispersion relation dE/dp do not differ very much from the free 
one. As a consequence the depth of the potential does not affect the spectra very much. 
Also possible minima of the dispersion relation |L4|, [15] do not affect the outside spectrum, 
because they also occur at energies smaller than the rest mass. 

In conclusion, we find, that in the (not very realistic) case of a static potential well, 
only small effects of the potential on the particle spectrum can be observed, too small to 
account for the observed enhancement. The main reason is that the particles which may 
escape the potential feel only a rather weak potential with very moderate momentum 
dependence. 

Of course as already pointed out, in reality the potential is not static because its 
source, the fireball, expands. Consequently the potential decreases as time continues so 
that eventually all particles may escape. The observed enhancement could, therefore, still 
originate from those particles, which would be trapped inside the static potential. 

This would be the case if these particles leave the system as soon as energy conservation 
allows the to do so, i.e once their energy becomes larger than their rest mass. However, 
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such a scenario, depends very much on the time scales involved in the problem. It would 
require that the expansion of the fireball is very slow compared to the velocity of the low 
momentum particles. Since the fireball is made out of pions as well this seems to be very 
unlikely. 

On the other hand, if the fireball expands with a velocity faster than the low mo- 
mentum particles these particles will essentially remain inside the potential well (fireball) 
until the potential has vanished. In this case only the fast particles would feel an effect 
of the potential as they have to climb a potential well of finite depth. At the soft part of 
the spectrum, however, we would not expect any effects of the potential^. 

The question of the time scales involved can best be answered in a model calculation. 
In the next section we, therefore, will study the full expansion of the fireball and the effect 
of the collective potentials on the spectra in a transport model. This approach should 
provide a reasonable simulation of the expansion and the time scales involved. 



4 This later scenario has been studied in great detail in case of anti-protons in ref. JlS] 
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4. Expansion of the fireball 



In order to study the dynamic effects of the collective potential we have extended a 
cascade model |19|| to include the propagation of the particles in a mean field. Follow- 



ing the standard procedure [16] the particle coordinates and momenta are propagated 



according to Newton's equations of motion 

dr p dU 
dt E dp 

= -<E (15) 
dt dr K ' 

where the mean field potential U is given by eqn. (|6]). The temperature which is needed 
in order to determine the mean field potential is calculated from the density of pions 
assuming local thermal equilibrium. The pion-density is extracted from the actual particle 
distribution. At every timestep the radial dependence of the density is fitted with a 
function of the form 

p( r± ) = A(l + b 2 r 2 ± ) exp(-6 2 r| ) (16) 

The parameters A and b are determined by the root-mean-square radius of the distribution 
and the total number of pions. In the longitudinal direction, on the other hand, for a 
given timestep the density is assumed to be constant between z min and z max , where z min 
and z max denote maximum distance in positive and negative direction where particles 
have materialized. Comparing with the actual density distributions these assumption are 
very well justified for the conditions we are dealing with and which we will discuss belowQ. 
The error introduced is certainly not larger than the one one would have when using a 
spatial grid and having to deal with large density fluctuations ||20|| . 

The cascade includes n, r], p and uj mesons but no baryons^. All mesons decay accord- 
ing to their empirical lifetime and decay channels. For the pion-pion scattering measured 
phase shifts are used while for all other elastic processes a constant cross section of 20 mb 
is assumed. The only inelastic process taken into account is 7r + 7r <-» p. In addition the uj 
is allowed to decay into three pions. While this may not take into account all the details 
of the meson-meson scattering it certainly provides enough accuracy in order to lead to 
a reasonably realistic expansion scheme. Finally the model does not respect the Bose 
nature of the mesons, i.e. we do not have any Bose enhancement factors incorporated 



in the collision integral as done e.g. in refs. [|T|, [13| . Here, we are rather interested in 
dynamical effects due to collective potentials. 

5 This is essentially a result of the Bjorken initial conditions we have imposed here (see below). 
6 At CERN energies the ratio of protons to n~ is about 1:6 for S + S collisions M). 
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We shall specifically be concerned with the central 200GeV/a 16 + Au data of 
the NA35 collaboration [|l|. This experiment measures negative particles and does not 
identify the pions explicitly. Usually one assumes a ~ 10% admixture of kaons, electrons 
and anti-protons. For simplicity, however, we neglect this fact and assume all negatives 
to be pions. 

The number of initial particles is determined such that the measured rapidity dis- 
tribution of pions is reproduced. We further impose Bjorken-like [^] initial conditions: 
rapidity yboost and longitudinal coordinate of a locally thermal distribution are uniquely 
related by 

z = r sinhy boost (17) 

where To may be considered a formation or hadronization time. In other words particles 
materialize, i.e. participate in the expansion, only after their proper time is larger than the 
hadronization time r . Because of eqn. ( |T7| ) this implies that in the cm. system particles 
which are created at a large longitudinal distance from the center will materialize at a 
later time. The initial radial distribution is assumed to follow the density profile of the 
oxygen projectile. 

We will study different ways of populating the initial momentum space. One possibility 
is to assume local thermal and chemical equilibrium so that with given density the initial 
local temperature and the admixture of higher resonances is fixed. Thus, in this approach 
the free parameters of the model are the initial number of particles and tq. Those can 
be determined by fitting the measured rapidity distribution and slope of the high energy 
part of the pion spectrum. The lowest hadronization time which still leads to acceptable 
transverse momentum spectra and rapidity distributions is r ~ 8fm/c. The number of 
initial particles is 380. 

The other procedure assumes that the momenta of all initial particles are distributed 
according to a slope parameter for p-p collisions of To = 135 MeV. Their relative multi- 
plicities are distributed proportional to their degeneracy, i.e. 7r:77:p:o; = 3:l:9:3 
These statistical weights are very close to what one obtains in string fragmentation models 
23fl . In order to reproduce the measured rapidity distribution we start with 280 particles 



in the initial state. 

The resulting rapidity distributions of pions for both initialization schemes obtained 
after 30 fm/c are shown together with the experimental data of the NA35 collaboration 
in fig. ||. These results have been obtained with the pure cascade without potentials. 
However, since in our model the density and hence the potentials are assumed to be 
constant along the longitudinal direction, the rapidity distribution is not affected by the 
potentials. 

Let us now turn to the effect of the collective potentials. In fig. |^ and we show the 
resulting transverse momentum spectra using the different initial conditions described 
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above. Fig. ^| corresponds to the local equilibrium initialization while fig. |7| represents 
the string fragmentation picture. Different expansion schemes are studied: 



1. The particles do not interact but may decay (full line). 

2. Particles may collide but do not feel any mean field forces (short-dashed line). 

3. Particles collide and are subject to the mean field forces (long-dashed line and 
dashed-dotted line). 

We display only the results for parametersets 1 and 2 of table [3], because the results ob- 
tained with parameterset 3 are indistinguishable from the one obtained with parameterset 
2. 

First let us point out the difference between the expansion with and without collisions 
(full and dashed histogram) in the low p t spectra exhibited in fig |7|. Since the initial 
conditions do not correspond to local chemical equilibrium the collisions are still very 
effective and provide a significant enhancement a low transverse momentum. This result 



should also be seen in connection with the findings of ref. ||13|1 , where the effect of Bose 
statistics in the collision integral is studied. Considering the fact that our results have been 
obtained without Bose statistics, we conclude that a great part of the effect demonstrated 



in ref. [0 can be explained without any Bose statistics. This supports the finding 
pointed out in the aforementioned reference, namely that the effect of Bose statistics in 
the final state and in the matrix element cancel each other to a large extent. In case of 
the local-equilibrium initialization (figure []) the effect of the collision is very small. One 
of the reasons simply is that we we had to chose a very high hadronization time To, in 
order to reproduce the high momentum piece of the spectra. As a consequence the initial 
density of particles is comparatively low and thus the collisions are not so effective. 

For either initial conditions the mean field leads only to very little enhancement at low 
Pi. Also the stronger mean field parameterization does not provide more enhancement. 
Thus, the expansion does not provide the additional enhancement compared to the static 
limit. Quite to the contrary, the effects found in the static limit are reduced because the 
time-averaged depth of the potential well is smaller as a result of the expansion. And 
certainly, the enhancement due to the pions being trapped in the static potential, as 
discussed in the last section, does not show up. Obviously, the velocity of the fireball 
expansion is faster than the velocity of the low momentum pions. As already mentioned 
this is actually what one should expect, because the relevant part of the fireball, which 
provides the potential, consists of pions. The average velocity of these pions, however, 
is much larger then the one of the soft pions. Therefore, as long as there are no exotic 
transport phenomena, we expect the expansion velocity of the source of the potential to 
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be higher than the velocity of the low momentum pionsQ. As a consequence the momen- 
tum of the pion remains essentially unchanged during the expansion. Thus, as long as 
we consider only one fluid (in our case the pions) an attractive potential does not lead 
to considerable enhancement in the spectrum at low momenta, because the velocity of 
the particles and the expansion velocity of the potential well are intimately related. This 
would be different if one had two different fluids, one which provides the potential, and 
one, the spectrum of which we want to study. In this case one could very well imagine that 
the expansion proceeds essentially adiabatically so that all bound pions would leave the 
potential with minimal momentum. Such a possibility actually exist at BEVALAC/SIS 
energy (~ 1 GeV) heavy ion collisions. There the pions feel a collective potential through 
the coupling to the delta-hole channel [p4| , which looks very similar to the one we have 
studied here. In this case, however, the potential is provided by the nucleon fireball. Since 
the mass of the nucleon is considerably larger than the one of the pions, it could very 
well be that the expansion velocity of the nucleon fireball is smaller than the velocity of 
the low momentum pions. Therefore, at BEVALAC/SIS energies the pion spectra could 
reveal information about the long sought in medium pion dispersion relation^. 



Conclusions 



In conclusion, we could show that the observed enhancement at low transverse mo- 
mentum in the spectra of pion cannot be accounted for by collective mean fields. We 
have pointed out, that this result depends very much on the time scales involved in the 
expansion. As long as the potential is provided by the same kind of particles as the ones 
the spectra are being studied of, we do not expect an enhancement, because in general the 
expansion velocity is faster as the velocity of the soft particles. This would be different at 
BEVALAC/SIS energies where the source for the pion collective potential would be given 
by the nucleons, which, as a result of their larger mass, would most likely expand slower 
than the soft pions. 

In the case where the particles initially have been distributed according to the mo- 
mentum spectrum of proton-proton collisions, we found that the inclusion of particle 
collisions leads to considerable enhancement at low transverse momentum in comparison 

7 It is actually very difficult to imagine that an attractive two-particle interaction should slow down 
slow particles even more. Let us consider a pair of pions, fast and slow, which interact by an attractive 
potential. Clearly the effect of the potential is, to slow down the fast pion and to accelerate the slow one, 
such that both at the end have the same velocity. 

8 Indeed pionic spectra at these energies exhibit an enhancement a low energies similar to those dis- 
cussed. This possibility will be addressed in a separate publication E5| . 
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to the free decay, which would more or less correspond to the simple folding of the p-p 
data. This effect, however does not fully account for the enhancement observed in the 
data. 



Taking into account the results of ref. [13] it seems that the soft pion puzzle seems still 



unresolved. The data require either a very short hadronization time (r — 1 fm/c) or the 
absence of mesons heavier than the pion. Both assumption certainly require a nontrivial 
scenario in order to be acceptable. 
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Figure captions 



Figure 1: Mean field potential for M res = 770 MeV for different temperatures (T = 
150, 200, 250 MeV) (full lines) together with fits according to eqn. (||) (dashed lines). 



Figure 2: Mean field potential @ for different resonance masses. Full line: M r 
500 MeV, Long-dashed line: M res = 770 MeV Short-dashed line: M res = 1100 MeV. 



Figure 3: Spectra for static potential based on eqn @ and parameterizations of table [l|. 
Free: full line; Set 1: short-dashed line; Set 2: long-dashed line; Set 3: dashed-dotted 
line. 

/ 
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Figure 4: Dispersion relation for parameterizations of table p]. Labels as in figure |3| 



Figure 5: Rapidity distribution for the two initialization schemes described in the text. 
Full histogram: local thermal equilibrium; dashed histogram: p-p distribution. 



Figure 6: Transverse momentum spectra for local equilibrium initialization: Full line: no 
mean-field, no collision; short-dashed line: no mean-field, collisions; long-dashed line: 
mean field (Set 1) collisions; dashed-dotted line: mean field (Set 2) collisions. 



Figure 7: Same as figure [|, but with p-p distribution in the initial state. 
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